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Abstract 

Background and Purpose: Recent data suggest that early symptoms may be related to cortex alterations in CADASIL 
(Cerebral Autosomal-Dominant Arteriopathy with Subcortical Infarcts and Leukoencephalopathy), a monogenic model of 
cerebral small vessel disease (SVD). The aim of this study was to investigate cortical alterations using both high-resolution 
T2* acquisitions obtained with 7 TesIa MRI and structural Tl images with 3 TesIa MRI in CADASIL patients with no or only 
mild symptomatology (modified Rankin's scale <1 and Mini Mental State Examination (MMSE) >24). 

Methods: Complete reconstructions of the cortex using 7 TesIa T2* acquisitions with 0.7 mm isotropic resolution were 
obtained in 1 1 patients (52.1 ±13.2 years, 36% male) and 24 controls (54.8±11.0 years, 42% male). Seven TesIa T2* within 
the cortex and cortical thickness and morphology obtained from 3 TesIa images were compared between CADASIL and 
control subjects using general linear models. 

Results: MMSE, brain volume, cortical thickness and global sulcal morphology did not differ between groups. By contrast, 
T2* measured by 7 TesIa MRI was significantly increased in frontal, parietal, occipital and cingulate cortices in patients after 
correction for multiple testing. These changes were not related to white matter lesions, lacunes or microhemorrhages in 
patients having no brain atrophy compared to controls. 

Conclusions: Seven TesIa MRI, by contrast to state of the art post-processing of 3 TesIa acquisitions, shows diffuse T2* 
alterations within the cortical mantle in CADASIL whose origin remains to be determined. 
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Introduction 

Cerebral autosomal-dominant arteriopathy witli subcortical 
infarcts and leukoencephalopathy (CADASIL) is a rare hereditary 
cerebral small vessel disease (SVD) caused by mutations of the 
N0TCH3 gene [1]. It is now widely recognized as a unique model 
for the study of more prevalent forms of SVD related to aging and 
hypertension [2] . The pathophysiology of SVD is thought to be 
mainly of subcortical origin, hence the term "subcortical ischemic 
vascular dementia" [3], but during the last decade the cerebral 
cortex was shown to be a key player in both sporadic forms of 
SVD and in CADASIL [4-8]. 

In a recent longitudinal study in CADASIL, processing speed 
slowing, but not global cognitive worsening nor increased disability 
was related to alteration of sulcal morphology, suggesting that 
early cognitive changes may be more specifically related to sulcal 



morphology than to other MRI markers [9]. The observed 
modifications of sulcal morphology could however be linked to 
underlying white matter structure changes [10] or to actual cortex 
involvement. The latter hypothesis could involve focal intra- 
cortical small-sized lesions [4], cortical demyelination [11] or 
increased iron content due to secondary degeneration, as reported 
previously in deep nuclei in CADASIL [12], aU of which may 
translate in alterations of T2* measured using 7 Tesla MRI within 
the cortical mantle [13]. 

The aim of the present study was to investigate the cortex 
morphology using state of the art post-processing of 3 Tesla 3D Tl 
acquisitions together with its structure using high-resolution T2* 7 
Tesla acquisitions in CADASIL patients at the early stage of the 
disease, with no or only mild symptomatology (modified Rankin's 
scale Si and MMSE S24), by comparison to age and sex 
matched controls. 
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Methods 

Participants 

The ethics committee of DRCD (Departement de la Recherche 
Clinique et du Developpement) of AP-HP (Assistance Publique- 
Hopitaux de Paris) validated the protocol and all subjects gave 
their written consent for participating in the study. The present 
data originate from a prospective study of CADASIL patients with 
genetically confirmed diagnosis, without dementia (MMSE score 
S24) and without significant disability (modified Rankin's scale £ 
1) dedicated to the study of early cognitive symptoms in this 
disorder. Twenty-three patients and 29 age and sex matched 
controls were included in the study. Sixteen patients and 25 
healthy controls had high quality 7 Tesla MRI covering the whole 
hemispheres and 3 Tesla MRI and were included in the present 
study. Clinical and demographic data were collected, including 
age, sex, smoking habits, alcohol intake, and body mass index. 
Laboratory evaluation (which included complete blood count, 
fasting glucose, high-density Kpoprotein, low-density lipoprotein, 
and total cholesterol levels) was performed in all subjects. AU 
participants underwent a comprehensive neuropsychological 
evaluation preceded by a clinical psychological interview. Global 
cognitive performances were assed using the Mini Mental State 
Evaluation (MMSE) [14]. Disability was assessed using the 
modified Rankin's scale (mRS, from 0- no disability to 5— major 
disability), in which scores inferior or equal to one denote no 
disability [15]. An executive clock drawing task, CLOX [16] was 
used as a planning task, with both spontaneous and copy subtests. 
A verbal fluency task (phonological and semantic) [17] was 
administered and working memory was assessed using the working 
memory index from the Wechsler Memory Scale Third Edition 
(MEM-III) [18]. Specific computerized test were performed with 
various levels of executive demand but are not detailed here as 
they are not relevant to the present study. Presence of depression 
was defined foUowing DSM-IV-TR criteria [19] and a depression 
rating scale, MADRS, was obtained for all subjects [20] . Subjects 
were considered apathetic if they had a Starkstein's Apathy Scale 
score &14 and fulfilled the Starkstein Structured Clinical 
Interview for Apathy criteria [21]. 

MRI protocol 

Subjects underwent the same day both 3 Tesla and 7 Tesla 
MRI acquisitions at NeuroSpin (CEA, Gif-sur-Yvette, France). Tp 
weighted 3D images were acquired with a 3 Tesla Tim-Trio MRI 

scanner (Siemens Healthcare, Erlangen, Germany) equipped with 
a 12-channel head coil, using a standard sagittal MPRAGE 
sequence (in plane resolution: 1 x 1 mm^, slice thickness =1.1 mm, 
TR = 2300 ms, TE = 2.98 ms, TI = 900 ms, FA = 9°, 
BW = 238 Hz/pixel, time of acquisition = 7 '45 min). Ultra-high 
resolution MRI acquisitions were performed on a 7 Tesla MRI 
scanner (Siemens Healthcare, Erlangen, Germany) equipped with 
a lTx/8Rx head coil (Rapid Biomedical, Wurzburg, Germany). 
2D axial gradient-echo sequences were acquired following AC-PC 
orientation, using seven or eight blocks depending on subject' 
brain size. The choice of 2D blocks rather than whole brain 3D 
acquisition was made because of the impossibility for most subjects 
to remain free of movements during the necessary- long acquisition 
periods. Each block comprised 20 slices with two different echo 
times using the following MR parameters: in plane resolution: 
0.7x0.7 mm^, sKce thickness = 0.7 mm, TR = 900 ms, 
TEi = 13.7 ms, TEg = 29.9 ms, FA = 65°, BW= 70 Hz/pixel, 
NA = 2, time of acquisition = 6 mn 46 s). Blocks were acquired 
parallel to each other with no overlap and were repeated, if 
possible, each time image was of poor quality. Total scan time by 



subject was about 1 hour on the 7 Tesla scanner. As FLAIR and 
T2''' images were already obtained within 6 months on a 1 .5 Tesla 
scanner (GE Healthcare, Milwaukee, Wisconsin) as part of an on- 
going prospective study in patients, theses sequences were used for 
the determination of white matter lesion volume and of number of 
microhemorrhages. 

Image processing 

Quantification of subcortical MRI markers. Identi- 
fication and quantification of MRI markers of SVD were 
performed according to recent standards [22]. White matter 
hyperintensity masks were automatically determined on all axial 
FLAIR slices from the base of the cerebellum to the vertex, using 
adaptative intensity thresholding [23]. Lacunes and dilated 
perivascular spaces were manually identified by experienced 
readers and lacune masks for each patient were automatically 
determined using simple connectedness algorithm with manual 
corrections when necessary [24]. All masks were manually 
validated in a second step by double reading, as previously 
reported [9]. Volumes were obtained by multiplying the number 
of voxels covered by the masks with their size depending on the 
sequence. The differentiation of lacunes from dilated perivascular 
spaces was based on current recommendations [22]. The number 
of cerebral microhemorrhages, defined as rounded foci smaller 
than 1 0 mm in diameter, hypointense on gradient-echo sequences, 
and distinct from vascular flow voids, leptomeningeal hemosider- 
osis, or non-hemorrhagic subcortical mineralization, was recorded. 

Quantification of global and sulcal morphology. AU 
images were processed using the Brain Visa suite (http:/ /www. 
brainvisa.info). Brain parenchymal fraction (BPF) was defined as 
the ratio of brain volume to intracranial cavity determined on the 
3D-T1 sequence, as previously reported in various disorders and 
in CADASIL [25]. Brain Visa was also used to reconstruct the 
cortical sulci [24]. Identification of sulci was made using a 
dedicated algorithm. For all subjects, results were carefully 
checked for both 2-dimensional segmentation masks and 3- 
dimensional reconstruction of sulcal shapes. Sulcal depth and 
sulcal width were obtained for 57 sulci on each hemisphere, and 
then averaged to determine for each subject mean sulcal width and 
mean sulcal depth. 

Cortical thickness analysis. Three Tesla Ti images were 
processed using FreeSurfer version 5.1.0 (http://surfer.nmr.mgh. 
harvard.edu), which produced surface-based data for each subject 
[26]. FreeSurfer pipeline includes data resampUng, standardized 
stereotaxic space normalization, intensity non-uniformity correc- 
tion using N3, skull stripping, tessellation of the gray/ white matter 
boundary, automated topology correction, and surface deforma- 
tion to optimally delineate the gray/ white matter and gray 
matter/cerebrospinal fluid boundaries [27,28]. As previously 
reported, masks of white matter hyperintensities were obtained 
from FLAIR sequences and registered to 3D Tl images [9,29]. 
Reconstructed datasets were always visually inspected for accura- 
cy. 

Cortical thickness maps were estimated by calculating the 
distance between the pial and white matter surface at each vertex 
from the 3 Tesla Tl -weighted image [28]. For group analysis, 
each surface was spherically inflated and automatically registered 
to a canonical spherical tf;mplat(; using cortical features. Thickness 
data from each subject were mapped to the average surface and 
smoothed using a Gaussian filter with a full-width half maximum 
of 10 mm before running general linear models (GEM). We 
contrasted CADASIL patients versus control subjects, covarying 
for age, gender and level of education. Thickness GEM results 
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were corrected for multiple comparisons using a false discovery 
rate (FDR) of 0.05 [30]. 

Post-processing of 7 Tesla acquisitions. As high-resolu- 
tion acquisitions at 7 Tesla were performed using multiple 
acquisitions of 2D slabs, a specific pipeline was set up in order 
to reconstruct the high-field whole brain image from those 
different blocks, as previously reported by another team in the 
context of multiple sclerosis [11]. The goal was to register the 
multiple 7 Tesla T2*-weighted blocks onto the 3 Tesla Tl- 
weigthed anatomical scans (see Fig. 1). 

a) Intra-subject registration. AH 7 Tesla T2* blocks were first 
concatenated to create a rough whole brain reconstruction (Fig. 1). 
We then used the FSL FLIRT algorithm (http://www.fmrib.ox. 
ac.uk/fsl) to perform an initial aSine registration between whole 
T2*-weighted volume at second echo time and 3D Ti volume 
which was over-sampled to an isotropic resolution of 0.7 mm. 
Next, boundary-based registration (BBR) was initialized with this 
latter transformation and ran separately for each block using the 
second echo time image for a better tissue contrast. BBR principle 
is to align the input image on the reference image by maximizing 
the intensity gradient across tissue boundaries [31]. BBR was 
particularly designed for aligning partial-brain images to whole- 
brain and for registering images with different contrasts. Quality of 
the block registration was visually checked systematically. If not 
satisfactory, manual adjustment was iteratively performed to 
change the initial conditions before BBR registration. Fine 
registered blocks were finally combined to create the reconstructed 
T2* volume. Data were dropped if registration still failed or if 
FreeSurfer pipeline was unable to reconstruct cortical surfaces. 

b) 7 Tesla volume pre-processing. Potential nuU voxels 
corresponding to spaces between blocks were excluded from 
analysis. The reconstructed 7 Tesla T2* volume at second echo 
time was normalized by the first echo time volume and then 
corrected for non-uniformity intensity using N3 algorithm [27]. 



This was achieved to limit large spatial variations of signal due to 
BO and Bl heterogeneities. 

c) 7 Tesla cortical signal analysis. Signal at the midline of the 
cortical ribbon as defined by the internal and external surfaces 
computed by FreeSurfer was extracted from the normalized 
volume and mapped onto each subject's cerebral surface. Each 
surface was then registered to the same average surface template of 
FreeSurfer with a 2D Gaussian smoothing (FWHM= 10 mm). 
Vertex-wise GLM was then run, contrasting CADASIL patients 
versus control subjects controlling for age, gender and level of 
education. Surface GLM results were corrected for multiple 
comparisons using a false discovery rate (FDR) of 0.05. Regions 
with significant differences were identified using an automatic 
parcellation of cortical gyri and sulci [32]. For each identified 
region, mean and standard deviation of 7 Tesla T2* signal and 3 
Tesla cortical thickness were further calculated in the native space 
of each subject. Global T2* signal within the cortical mande was 
also calculated for each subject in its native space. High quality 
data was obtained for 1 1 patients and 24 controls. 

Statistics 

Statistical analyses reported in Table 1 were made using the R 
software (http://www.r-project.org/). For categorical variables, 
Chi2 tests were used. For numerical variables, Wilcox-tests were 
used given the small sample sizes. 

Results 

Characteristics of both groups included in the analysis are 
detailed in Table 1. The two groups did not differ in terms of age, 
gender or MMSE (29 for median MMSE in both groups). BPF, 
mean sulcal width and mean sulcal depth did not either differ 
between groups. Patients difiered relatively to controls in the 
working memory task and tended to have higher scores on apathy 
and depression scales. 




Intra-subject registration 




Figure 1. Image processing pipeline for surface-based intra-subject and inter-subject registration. 1. Original 7 Tesla T2*-weighted 
blocks (axial image of a CADASIL patient is shown for example) are concatenated. 2. Reconstructed volume is registered to the 3 Tesla T,-weighted 
image to produce an initial alignment. 3. Each 7 Tesla block is then registered to the 3 Tesla image using the boundary-based registration algorithm 
initialized with the latter transformation. All blocks are combined to get a whole-brain 7 Tesla T2* reconstructed volume. 4. Pre-processing include 
signal normalization by the first echo time image intensity, correction for intensity inhomogeneity and exclusion of null points corresponding to 
voids between blocks. 5. 7 Tesla normalized T2* signal is mapped onto each individual's cortical surface at the midline of the cortical ribbon to study 
cortical signal alterations. 6. Each surface is smoothed and registered to a common template. 
doi:1 0.1 371/journal.pone.01 0631 l.gOOl 
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Table 1. Characteristics of CADASIL patients and control subjects. 









CADASIL patients (N = 1 1 ) 


Healthy controls (N = 24) 


P- 

value 


Age, mean ± sd, range 




52.1 ±13.2, 32.1-70.6 


54.8±11.0, 30.1 


-71.4 


0.57 


Male gender, n (%) 




4/11 (36%) 


10/24 (42%) 




0.99 


Right handed, n (%) 




11/11 (100%) 


24/24 (100%) 




0.99 


Level of education (years) 




11.4±3.5 


13.6±3.6 




0.10 


Neuropsychological tests 


MMSE, mean, median, range 




28.5, 29, 24-30 


29.0, 29, 26-30 




0.45 


Fluency letter p 




24.1 ±9.8 


25.2±7.9 




0.82 


Fluency category 




22.1 ±6.1 


22 .3 ±4.4 




0.75 


Clox copy 




14.8 ±0.4 


14.8 ±0.4 




0.84 


Clox spontaneous 




12.8±1.5 


13.3±1.4 




0.18 


Working memory index 




89.8 ±11.4 


101.3±10.0 




0.005 


Apathy scale 




9.6±3.1 


6.8±2.9 




0.03 


Depression score 




4.9±5.7 


1.5±2.9 




0.02 


Clinical data 


SEP (mmHg) 




122.3 ±10.5 


128±14.8 




0.35 


DBP (mmHg) 




75.6±5.8 


77.4±14.5 




0.90 


Active smoking, n (%) 




1/11 (9.1%) 


4/24 (16.7%) 




0.94 


Glucose levels (g/l) 




1.03±0.13 


0.95 ±0.07 




0.07 


LDL-cholesterol (g/l) 




1.09±0.29 


1.27 ±0.29 




0.18 


HDL-cholesterol (g/l) 




0.62±0.13 


0.64±0.18 




0.88 


Triglycerides (g/l) 




1.30±1.04 


1.03 ±0.68 




0.57 


Imaging markers 


EPF 




0.85 ±0.03 


0.84±0.03 




0.32 


Mean sulcal width, mean ± 


sd (mm) 


1.36±0.20 


1.45 ±0.24 




0.23 


Mean sulcal depth, mean ± 


sd (mm) 


13.41 ±0.84 


13.76±0.47 




0.22 


White matter lesion volume. 


mean, median, range in mm^ 


65511.6, 47544, 7249.8-233797.3 


No control had 


significant lesions 


NA 


Volume of lacunes, mean, median, range in mm^ (n = 5/11, 45%)^ 


423.5, 113.1, 14.0-1174.2 


No control had lacunes 


NA 


Number of microhemorrhages, mean, median, range (n = 3/1 1, 27%)^ 


2.3, 1, 1-5 


No control had 


microhemorrhage 


NA 



MMSE: Mini Mental State Examination, BPF: Brain Parenchymal Fraction, 
NA: not applicable. 

doi:l 0.1 371 /journal.pone.Ol 0631 1 .tOOl 



'in patients with such lesions (number given in parentheses). 



7 Tesia cortical nneasures 

Mean normalized T2* signal for both groups is represented in 
Fig. 2. These cortical maps highlight regional differences across 
the cortex. While sensorimotor, visual and auditory areas 
exhibited lower values, higher values were found in the middle 
and superior frontal areas and in the cingulate cortex. Large 
hypointense inferior regions in frontal and temporal cortices were 
due to susceptibility effects. This pattern of T2* changes was found 
in both groups but differences between regions were visually 
smaller in CADASIL patients who showed a general increase in 
T2* in all cortical areas compared to control subjects. Mean T2* 
signal in the whole cortex was also found significantly higher in 
patients than in controls (0.73±0.07 vs 0.68±0.13, p = 0.01 after 
adjustment for age, gender and level of education). 

GLM results contrasting both groups showed several regions 
with a significant T2* increase in CADASIL patients compared to 
controls after adjustment for age, gender and level of education. As 
represented on the average pial surface, differences were largely 
symmetrical (Fig. 3). Table 2 reports significant local T2* changes 
between groups. There was no region where controls showed 
increased T2* compared to CADASIL patients. Higher T2* in 



patients with CADASIL were essentially located in parietal, frontal 
and occipital cortices. The left hemisphere showed more 
significant differences between groups. No significant correlation 
between T2* increase and BPF, white matter lesion volume, 
volume of lacunes or number of microhemorrhage were detected 
(data not shown). Comparing the 1 1 CADASIL patients with 1 1 
age and sex matched controls yielded similar results with only 
reduced significance (see Figure SI). Mean cortical thickness in 
regions showing a difference in 7 Tesla T2* were not statistically 
different (Table 2). 

Cortical thickness 

GLM results contrasting both groups showed no significant 
differences between CADASIL patients compared to controls after 
adjustment for age, gender and level of education, except a single 
area showing reduced thickness in patients in the precentral gyrus 
of the right hemisphere (Fig. 3). Mean thickness over the whole 
cortex was not different between controls and CADASIL patients 
(2.59±0.14 mm vs 2.57±0.09 mm, p = 0.49 after adjustment for 
age, gender and level of education). 
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Figure 2. Normalized T2* signal in CADASIL and control subjects. Normalized T2* signal assessed at 7 Tesia across both groups is displayed 
on an inflated average surface for the right hemisphere. In both groups, hypointense signal is found around the central sulcus (1) due to high myelin 
content. As also previously reported, control subjects exhibit shorter T2* in visual (2) and auditory (3) areas and longer T2* in superior and middle 
frontal sulci (4) and in the posterior (5) and anterior (6) cingulate. Temporal (A) and frontal (B) signal loss due to susceptibility effects. CADASIL 
patients show an overall longer T2* compared to controls, leading to less variations between different cortical regions. 
doi:1 0.1 371 /journal.pone.01 0631 1 .g002 



Relationships between T2* and clinical status 

None of the neuropsychological variable was found significantly 
associated with T2* cortical signal using FDR threshold. Some 
marginal associations were observed using a more permissive 
threshold (p<0.05) as represented in Figure S2. 

Discussion 

In the present study, we observed a strong and diffuse T2* 
increase within the cortical mantle in CADASIL patients with no 
or only mild symptomatology by comparison to age and sex 
matched controls. By contrast, conventional state-of-the art post- 
processing of high-resolution 3D T 1 images showed only minor, if 
any, group differences. The increased T2* was significant in 
various frontal and parietal regions most often following a 
symmetrical pattern on both hemispheres. These T2* alterations 
were not related to cortical thinning. 

The observed increase of T2* at ultra-high field in patients 
could have several explanations. Indeed, T2* signal mostly relies 
on proton density and spin phase coherence between protons. 
When proton density or spin phase coherence increase, then the 
T2* signal is increased. The presence of intracortical infarcts, as 
previously depicted pathologically in a CADASIL patient [4], 
could lead to an increased T2* through an increase of proton 
density because of the presence of intracortical fluid-filled cavities, 
but this hypothesis has been ruled out by visual inspection 
although we cannot formally exclude the presence of numerous 
cortical microinfarcts not visible on MRI. However, this appears 
very unlikely given that the patients were all at the early stage of 
the disease. 

The increase of T2* could also be secondary to a reduction of 
the loss of spin phase coherence between protons. In the cortex, 
both myelin and iron contribute to the T2* relaxation, with a 



predominant role of iron [33,34]. Decreased iron concentration or 
demyeUnation could be involved in the T2* increase according to 
MR tissue contrast modeling from both histology and post- 
mortem MRI scanning [34] . By contrast, secondary degeneration 
related to the accumulation of lacunar lesions in subcortical areas, 
as reported on conventional imaging recently [10], is unlikely 
explaining our results as it would rather lead to a decrease of the 
T2* signal because of an increased iron concentration as 
confirmed pathologically in basal ganglia [12]. 

As deoxygenated red blood cells are paramagnetic and induce 
loss of spin phase coherence, another possible origin of these T2* 
increase within the cortex might be the variations in venous 
vascular density and/or blood oxygenation [35]. A reduction of 
vascular density, in line with observations reported in a mouse 
model of CADASIL would imply less spin dephasing and thus 
increased T2* [36]. Capillary loss has also been observed 
pathologically in the cortex of patients with extensive white- 
matter hyperintensities [37], and we have recentiy reported using 
7T MRI an increase of T2* within the white matter in CADASIL 
patients having less visible veins in the centrum semiovale [38] . 

Otherwise, myelin-associated iron may also be a major source of 
intracortical magnetic susceptibility-based contrast [33]. In mul- 
tiple sclerosis, T2* increase observed at 7 Tesla have been 
associated with general demyeUnation [11]. Edema may also be 
involved in the observed T2* increase. We have reported that 
brain volume was increased in CADASIL patients with high loads 
of white matter hyperintensities [39] and a recent study in a mouse 
model of CADASIL have shown an intramyelinic edema as soon 
as the early stage of the disease, without blood-brain-barrier 
involvement [40]. Other potential sources of increased T2* at 
ultra-high field include inflammation and gliosis [41], but to our 
knowledge, they have not been reported to be diffuse in the cortex 
of CADASIL patients. 
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3T cortical thickness decrease in Patients > Controls 
P-value 0.001 0.00001 



7T T2' increase in Patients > Controls 
P-value 0.003 0.0001 




1 : middle occipital gyrus and sulcus 

2 : intraparietal sulcus and 
transverse parietal sulci 

■ angular gyrus 
postcentral sulcus 
middle frontal gyrus and sulcus 
superior frontal gyrus 
paracentral gyrus and cingulate 
sulcus 

8 : precuneus j 



Figure 3. Surface-based GLM results: CADASIL vs Controls. GLM significance nnaps (p-value) overlaid on the average pial surface for both 
hemispheres on lateral (top row) and nnedial faces (middle row) for CADASIL patients (N = 1 1 ) versus control subjects (N = 24). 7 Tesia T2* and 3 Tesia 
cortical thickness differences were mapped in warm and cold colors, respectively. Results forT2* increase were also mapped for visualization on the 
right average inflated surface with cortical parcellation labeling (bottom row). There was no region in which CADASIL patients showed a significant 
reduced T2* compared to control subjects. By contrast, T2* was significantly higher in CADASIL patients in middle occipital gyrus and sulcus (1), 
intraparietal sulcus (2), angular gyrus (3), postcentral sulcus (4), middle frontal gyrus and sulcus (5), superior frontal gyrus (6), paracentral gyrus and 
cingulated sulcus (7) and precuneus (8). Cortical thickness was found reduced in patients only in the right precentral gyrus. Thresholds were set to 
FDR <0.05. 

doi:1 0.1 371 /journal.pone.01 0631 1 .g003 



In line with previous data, which showed alterations of the 
cortex morphology that parallel the clinical worsening in 
CADASIL [9], our findings emphasize the role of the cerebral 
cortex in CADASIL and more generally in SVD. Our results 
suggest that cortex alterations may appear early during the course 
of the disease, when patients stiU have no or only mild 
symptomatology and before global brain atrophy or cortical 
atrophy could be demonstrated using conventional imaging. By 
contrast to previous works, we did not observe cortical thinning in 
patients compared to controls [10], but in the present study, 
patients were less severe. Moreover, by comparison to previous 
studies in CADASIL, patients had lower loads of lacunar lesions in 
the present study, which may explain why cortical thinning could 
not be detected. The single area with decreased cortical thickness 
in right precentral gyrus should be interpreted cautiously given the 
multiple testing approaches in vertex wise analyses and this result 
may be due to chance alone. 

In this study, we did not observe any significant relationship 
between clinical variables and the increase of T2* signal in the 
cerebral cortex when taking into account multiple testing. When 
using a less conservative threshold, only minor associations were 
detected (Figure S2). Given the clinical diflFereiices between 
patients and controls, which were very small, it would have been 



surprising to find strong correlations with the clinical status. 
Further studies will be needed to determine whether these 
alterations are clinically meaningful. 

The main limitation of our study is the small number of patients 
as a major source of data loss was the incapacity to undergo long 
lasting acquisition (more than one hour in the 7 Tesia scanner). 
Importandy, results were left unchanged when the 1 1 patients 
were compared to 1 1 age and sex matched controls, with only a 
slight decrease in the level of significance (see Figure SI). Thus, our 
results could unlikely be explained by chance or selection biases. 
Moreover, the post-processing of images used in the present study 
is complex, involving multiple steps for both 3 Tesia and 7 Tesia 
acquisitions, which could have induced various biases and errors. 
Natives images were thus carefully checked and T2* measured at 
multiple positions, confirming the clear tendency to T2* increase 
in patients. Moreover, as patients did not show measurable global 
or cortical atrophy, these potential sources of systematic errors can 
be realistically excluded. 

Our study has also several strengths. High-field 7 Tesia T2* 
acquisition combined to an isotropic 0.7 mm'' resolution produced 
unique data to assess cortical T2* variations with reduced partial 
volume effects. The block acquisition protocol enabled to rescan 
individual blocks during the imaging session if subjects moved. 
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The surface-based intra- and inter-subject registration strategy 
made possible to overcome individual anatomical variability to 
probe global cortical signal changes. 

In summary, the use of 7 Tesla MRI high-resolution 
acquisitions allowed demonstrating strong and early involvement 
of the cortical T2* in CADASIL, undetected using state of the art 
post-processing analyses of conventional 3 Tesla MRI images. 
Further studies are needed to understand the underlying 
mechanisms of these alterations and their clinical consequences 
as well as to determine whether our findings can extend to more 
frequent forms of small vessel disease of the brain. 

Supporting Information 

Figure SI Effect of smaller control sample on surface-based 

GLM results. GLM significance maps (p-value) overlaid on the 
a\ erage left pial surface for the main analysis (right, 1 1 patients vs. 
24 controls) and for a matched number of controls (left, 1 1 patients 
vs. 11 controls). 
(TIFF) 
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